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In the note on closed orbit [1] , I mentioned that in the presence of imper- 
fections, an electron would follow a new path, called the closed orbit. Why 
can't it follow another path? Of course it can. If the electron receives a large 
force, it can even leave the accelerator. What I would like to show in this 
note is: 

If the electron is near the closed orbit, and the radiation loss is small, 
then the path of electron would get closer and closer to the closed orbit. 
To deomnstrate this, we start by noting that: 

1. If an electron happens to be on one point of the closed orbit - with the 
exact displacements and momenta - then it follows the closed orbit. 
This is because its path is determined by the same transfer matrices 
that give the closed orbit. 
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2. If an electron happens not to be on any point in the closed orbit, 
then it does not follow the closed orbit. It does not even intersect the 
closed orbit, otherwise the previous point applies and it would follow 
the closed orbit. 

We want to show that an electron not on the closed orbit would still 
approach the closed orbit after some time. In order to demonstrate this, 
the transfer matrices must first be linearised around the closed orbit. This 
means: 

1. finding a new, one turn transformation matrix for which the close orbit 
is now the origin 

2. this matrix would give the correct co-ordinates after one turn for small 
deviations from the closed orbit. 

Denote a point on the closed orbit at position s by (xo, x' Q , yo, y' 0 , zo, So, 1). 
Denote the deviation from this point of an electron at s by (x, x', y, y', z, 5, 1) 
With respect to original reference path, the actual co-ordinates of the 
electron is (x 0 + x, x' 0 + x', y 0 + y,y' 0 + y', z 0 + z, 5 0 + 5, 1). 

Denote the co-ordinates after one turn by (x 0 + x 2 , x' Q + x' 2 , y 0 + y 2 , y' 0 + 

y' 2 ,z 0 + z 2 ,S 0 + fa, 1). 

Start with the example of a horizontal bending magnet [5]. The 7x7 
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transfer matrix in the original reference path co-ordinates is given by 
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In terms of the closed orbit and the deviation, it can be written as: 
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For the point at s on the closed orbit itself, the initial and final positions are 
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the same, by definition of the closed orbit: 
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In order to linearise around the closed orbit, we subtract one matrix equation 
from the the other. To make it easier, we expand each matrix equation into 
the seven equations. As a starting example, we look at the second of the 
seven equations in eq. (2): 



x' 0 = x' 0 + m 26 S 0 + m 27 

and in eq. (3): 

x' 0 + x' 2 = (x' 0 + x) + m 2e (5o + S) + m 27 
Subtracting, we get 

x' 2 = x + m 26 5 



(4) 



(5) 



(6) 



Repeating this subtraction for all the equations in eqs. (2) and (3), we 
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get 
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Notice that the last row and last column of the matrix are all zero. So the 
matrix has contracted to a 6x6 matrix. The equation can be written in terms 



of the 6x6 matrix: 
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Note that the displacement vectors has only 6 elements now. 

In the same way, the linearised 6x6 matrices for all the other ring elements 



can be obtained. 



There is a slight complication for sextupole because the matrix elements 



depend on x and y [2]: 
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This can be expanded into 7 equations as before. The 2nd equation is: 

x[ = x' 0 + m 26 5o + m 27 (10) 

Substituting the expressions from [2] gives: 

x\ = x' 0 + l -\{xl + y 2 0 )5 0 - l -\{xl + yl) (11) 

In terms of the co-ordinates for the closed orbit and the deviation, this be- 
comes: 

x'q = x' 0 + ^X(x 2 0 - y 2 0 )5 0 - -X(x 2 0 - yl) (12) 

and 

x' Q + x' 2 = (x' 0 + x') + ^\[(xo + x) 2 -(yo + y) 2 }(S 0 + S)-^\[(x 0 + x) 2 -(y 0 + y) 2 } 

(13) 



subtract - neglect higher order such as xS, x 2 S, x 2 , yS, y 2 5, and y 2 , we get: 

x' 2 = x'+ ^X(x 2 0 - yl)8 - \x 0 x(l - So) + \y 0 y(l - S 0 ) (14) 

In [3], the factor of (1 — So) is apparently neglected, giving 

x' 2 = x' + ^X(xl - yl)5 - Xx Q x + Xy 0 y (15) 

Likewise, the 4th equation from eq. (9) is: 

y[ = y' 0 + m 46 5 0 + m 47 (16) 

Substituting the expressions from [2] gives: 

y[ = y' 0 - Xx 0 y 0 5 0 + \x 0 y 0 (17) 

In terms of the co-ordinates forthe closed orbit and the deviation, this be- 
comes: 

Vo = Vo - Xx 0 y 0 S 0 + Xx 0 y 0 (18) 

and 

Vo + V2 = Wo + y') - + x)(y 0 + y)(S 0 + S) + \(x 0 + x)(y 0 + y) (19) 

subtracting and neglecting higher order therms like x5, xy5, xy, and yS gives 

y' 2 =y' - \x 0 y 0 5 + \x 0 y(l + S 0 ) + \y 0 x(l + 5 0 ) (20) 

Again, in [3], the factor of (1 — So) is apparently neglected, giving 

y' 2 = y' - \x 0 y 0 5 + \x 0 y + \y 0 x (21) 
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Repeating the above steps for the 5th equation from eq. (9) gives: 

3 

z 2 = --\(xl-yl)x + 3\x 0 y 0 y + z (22) 

Note that in [3], the factor of 3 is not there. 

Combining the above results, the linearised sextupole transfer matrix is: 
/ 
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Note difference from expression in [3]: 
/ \ / 
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After obtaining all the 6x6 linearised matrices, the next step is to multiply 
them around the ring to get the one turn 6x6 transformation matrix at T{s). 



Then we solve the eigenvalue problem: 

T{s)E{s) = XE{S) (25) 

To see how this is related to the deviation from the closed orbit, write the 
eigenvalues as: 

A fc = exp(-a fc + ±i2nv k ) (26) 

for the eigenvectors E k (s) and E k (s)*, for k = I, II, III. The deviation from 
the closed orbit can be expressed in terms of the eigenvectors: 

X(s)=Y,CkE k (s) + c* k E k (sY (27) 
k 

After one turn, this would become 

X(s + 2nR) = ]T e~ ak [c k E k (s)e +l2 ™ k + c* k E k {s)*e~ i2 ™ k ] (28) 
k 

So if all three of the a k are positive, the deviation would get smaller and 
smaller with every turn around the ring, and the electron would approach 
the closed orbit. 

We now need to understand why or how a k could positive. 

It turns out that the damping depends on both the radiation loss at 
magnets, and energy gain in cavities. Consider energy gain u in a cavity. 
This causes the gradients x' and y' to decrease by a factor of 1 — u/E 0 . The 
reason is illustrated in fig. 1. 
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Figure 1: Damping effect of the cavity. 

Let Po = E 0 /c be the full electron momentum. Let P x be the actual 
transverse momentum. Then the gradient x' is given by 



(29) 



° ~ E 0 /c 

Figure 1 shows that if there is an energy gain of u, the gradient becomes 



x[ 



(30) 



E 0 /c + u/c 

Since the energy gain in the cavity u is small compared to the electron energy 
E 0 , a first order expansion gives: 



/ Px /, u 
1 E 0 /c V E 0 



1 - 



u 

En 



Xn 



(31) 



This transformation from x' 0 to x\ can be included in the 6x6 matrix by 
changing the value of 17122 from 1 to 



m 2 2 = 1 



u 
En 



(32) 



Applying the same reasoning to the gradient y' gives 



7B44 = 1 



u 
En 



(33) 
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Next, we consider the additional radiation loss in a magnet, as a result 
of the deviation 5 from the closed orbit. 

Previously, when the 7x7 matrix is used, the radiation loss is included 
in the 7th column. E.g. m 67 for the bending magnet [5] is used to include 
radiation loss in the value of Si. 

For the 6x6 matrix, this can be done by modifying the value of the diag- 
onal element m 66 for a magnet. The required change is sketched out in [4]. I 
shall fill in some details here. 

Following the approach of [4], the radiation loss is writtien in this form: 

„ 2 fe 4 B 2 \ ( E \ 2 

P=o (— »] (34) 



3 \m 2 c 3 J Kmc 2 , 

where m is the electron rest mass m e . This is derived from Lienard's formula 
[6] 

P = ^7V (35, 

where q is the electron charge e, a is the acceleration. Equation (34) derived 
from Lienard's formula using the following: 

E = •ymc 2 (energy mass relation) (36) 
r yma = Bee (Newton's second law and magnetic force) (37) 

We need to apply eq. (34) to a magnet element. In this element, the 
magnetic field B is fixed. We need to find the additional radiation loss when 
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energy E deviates from its value on the closed orbit by S For the relative 
amount S, the actual amount is ES. 

The radiation loss on the closed orbit at the magnet elment is e = Pt, 
where t is the time taken through the element, This amount has already 
be included in the 7x7 matrix when calculating the closed orbit When the 
energy is changed by ES from the closed orbit, there is additional radiation 
loss. We can find this in terms of e. 

Notice from eq. (34) that P is proportional to E 2 . If the energy changes 
by a small amount AE, P changes to P(l + 2AE/E). So when the energy 
change is ES, this become e(l + 2S). This gives a relative change of e(l + 
2S)/E 0 with respect to the reference energy E 0 . This energy loss has to be 
subtracted from 5. 

Split the energy loss into two parts: e/E 0 and 2eS/E 0 . The first part, 
e/E 0 , has already been subtracted during the closed orbit calculation. E.g. 
for the bending magnet, it is the term m e7 [5]. The other part, 2eS/E 0 , must 
now be subtracted from 5 for the deviation from closed orbit. 

That is, S must be changed to S — 2eS/E 0 . This can be done by modi- 
fying the element m§§ in the 6x6 transfer matrix for the magnet. Lets write 
S — 2eS/E 0 as (1 — 2e/E 0 )S. So when the transfer matrix is applied to the 
deviation vector, we need S to change to (1 — 2e/E 0 )S. We can achieve this 
if we replace the value of 1 for m 66 by (1 — 2e/E 0 ). 
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Summarising the results so far- Of the 6 diagonal elements in the 6x6 
transfer matrix, 3 needs to be modified: 

In a cavity, for the co-ordinates x' and y', 

m 22 = l-u/Eo (38) 
777.44 = 1 — u/Eq (39) 

In a magnet, for the co-ordinates 5, 

m 66 = 1 - 2e/E 0 (40) 

These expression give us a way to relate the energy gain and loss to the 
eigenvalues, and then to the damping rates. 

The product of the eigenvalues is equal to the determinant of the matrix. 
For an infinitesimal element of the ring, the determinant can be calculated 
using just the above elements. An infinitesimal element means a very short 
element in the ring. This can be a hypothetical fraction of a magnet, a 
fraction of a drift space, or a fraction of cavity. As explained in [4], "the 
only terms in the determinant which will be first order in the length of the 
element will be due to the diagonal terms of the matrix." 

To see that this is so, an inspection of the off diagonal elements for the 
transfer matrices given in Table 1 of [3] would show that they are indeed all 
to first order in length /. This is less obvious for the matrix elements for the 
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rf cavity, because it is given in terms of the potential difference V across the 
whole cavity. However, if we imagine dividing the cavity up into fractions, 
then V would also change according to the length of each fraction. 

A transfer matrix for an infinitesimal ring element, therefore, has the 
property that all the off diagonal matrix elements are first order in length. 
The diagonal matrix elements are 1 except 17122, rn^ and m 66 above. For 
these matrix elements, the difference from 1 comes from the energy gain or 
loss, which are also proportional to the length I of the ring element. 

To see that the determinant will depend to first order on the diagonal 
terms only, consider the following simple example: 



(41) 



1 + kid k 2 a 

A = 

^ k^a 1 + k 4 a 
where the terms in a are all much smaller than 1. The determinant is 



det(A) = (1 + ha)(l + k 4 a) - k 2 k 3 a 2 (42) 
= 1 + k±a + k^a + kik^a 2 — k 2 k 3 a 2 (43) 

To first order in a, the determinant is given by 1 + k±a + k^a. These include 
just the terms in the off diagonal elements, as pointed out in [4]. 

Likewise, for the 6x6 transfer matrix, the determinant would be given by 
1 + 2 8nm where 5 nn is the difference of the n th diagonal element from 1. For 



14 



a cavity, where 



u 

m 22 = 1 - (44) 

U , . 

m 44 = 1 --=- (45) 
the determinant would be 1 — 2u/E Q . For a magnet, where 

the determinant would be 1 — 2e/E 0 . 

Next, we find the determinant of the matrix for one revolution. This 
is obtained by multiplying all the transfer matrices of the elements in the 
ring. We use the property that the determinant of a product is equal to the 
product of the determinants: 

det(AB) = det(A) det(B) (47) 

Consider the product of a magnet transfer matrix and a cavity transfer ma- 
trix. Using the determinants obtained above, the determinant of the product 
is, to first order, 

(l 2M ) (l 2e ) ~ 1 2U 26 (48) 
V EqJ \ EqJ Eq Eq 

i.e. we just have to add 1 to all of the 5 nn in both matrices. This is readily 

extended to the product of all the matrices in the ring: 

det(T( S )) « 1 - ^ - ^ (49) 
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where X u is the total energy gain Uq around the ring, and X e is the total 
radiation loss. For the equilibrium condition, when the electron follows a 
steady orbit - the total energy gain is equal to the total radiation loss. Then 

det(T( S )) « 1 - ^ (50) 

To relate to the eigenvalues, we find the determinant from eq. (26): 

det(T(s)) = exp[-2(a 7 + a n + a in )\ (51) 

For small radiation loss, we expect damping to be small, so that 

det(T(s)) w 1 - 2(aj + a H + a in ) (52) 

Comparing eqs. (50) and (52), 

cui + an + oliii = (53) 

This is the Robinson's sum rule. 

In the limit of no closed orbit distortion and no coupling (the motion in 
one axis does not affect motion in another axis), the transfer matrice can be 
diagonalised into three 2x2 matrices, one for each axis. The above reasoning 
can be applied to each direction separately. The x and y directions depend 
only on magnet elements: 

ax = ay = w 0 ^ 
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The z direction depends only on rf cavities: 

a z = -ET~ (55) 

So for non zero radiation loss U 0 , the damping rate for all directions is posi- 
tive, and an oscillating electron would fall toward the reference path. 

However, even in the ideal case of no misalignment, there would be cou- 
pling. For instance, energy gain in a cavity would lead to decrease in 5, which 
would lead in turn to a change in x' through the next magnet. This means 
coupling between the x and z directions. 

However, for small radiation loss, and for small misalignments, we may 
expect that the values of the damping rates would not deviate too far from 
a x , a y and a z . Then the electron should approach the closed orbit after a 
number of turns. 

The Robinson's sum rule tells us that the sum of the damping rates 
is positive. It does not say whether a particular damping rate has to be 
positive. However, the damping rates can be calculated from the one turn 
transformation matrix. This provides a way to design the ring elements to 
make all three damping rates positive. The existence of operational rings 
where stable beams are achieved shows that this is possible. In such cases, 
the electron would indeed approach the closed orbit. 
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